A Model for Quantum Stochastic Absorption in Absorbing Disordered Media 
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Wave propagation in coherently absorbing disordered media is generally modeled by adding a 
complex part to the real part of the potential. In such a case, it is already understood that the 
complex potential plays a duel role; it acts as an absorber as well as a reflector due to the mismatch of 
the phase of the real and complex parts of the potential. Although this model gives expected results 
for weakly absorbing disordered media, it gives unphysical results for the strong absorption regime 
where it causes the system to behave like a perfect reflector. To overcome this issue, we develop a 
model here using stochastic absorption for the modeling of absorption by "fake" , or "side" , channels 
obviating the need for a complex potential. This model of stochastic absorption eliminates the 
reflection that is coupled with the absorption in the complex potential model and absorption is 
proportional to the magnitude of the absorbing parameter. Solving the statistics of the reflection 
coefBcient and its phase for both the models, we argue that stochastic absorption is a potentially 
better way of modeling absorbing disordered media. 
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I. INTRODUCTION 

Quantum and classical wave propagation and localiza- 
tion in disordered media is a well studied problem [1-3]. 
Recently in the literature, more attention has been paid 
to studying the effect of localization on coherent absorp- 
tion/amplification in coherently absorbing/ amplifying 
disordered media, both theoretically [4-11] and exper- 
imentally [12-18]. Light is Bosonic, so it can be am- 
plified or absorbed. On the other hand, electrons are 
Fermions, which can not be amplified; however, they can 
be absorbed in the sense of phase de-coherence. Co- 
herent absorption/amplification in a coherently absorb- 
ing/amplifying medium is a non-conserving scattering 
process where temporal phase coherence of a wave is pre- 
served despite absorption/amplification. Coherent back 
scattering (CBS) in a disordered medium, the main cause 
of weak and strong localization, is not affected by the ad- 
ditional presence of a coherently absorbing/ amplifying 
medium due to the persistence of phase coherence of the 
interfering waves. It was shown [4] that the localization 
can enhance the coherent light amplification in a coher- 
ently amplifying (i.e. lasing) disordered medium, result- 
ing in a self-sustaining mirror-less lasing action where 
coherent feed back for lasing action is supported by dis- 
order induced localization. In that same frame work, we 
also studied the effect of localization on the coherently 
absorbing disordered media and showed the similarities 
of the two processes for effective localization and effec- 
tive active lengths. This problem studied further in detail 
by different approaches [6-11]. The general approach to 
model coherent amplification and absorption is to add 
a constant complex part to the real part of the poten- 
tial. Depending on the sign of the complex potential, the 
model shows amplification or absorption. Modeling ab- 
sorption/amplification by a complex potential, however, 
always gives a reflection part due to the mismatch be- 
tween the real and the imaginary potentials. Studying a 



delta scatterer, which has both real and complex parts, it 
was shown [19] that absorption is not a monotonic func- 
tion of the strength of the complex potential. Although 
complex potential model works well for the case of weakly 
absorbing/amplifying disordered media, it gives unphys- 
ical results in the strong absorption/amplification regime 
because the system behaves like a perfect refiector. To 
address this issue, it was derived and demonstrated [5] 
that absorption without reflection can be modeled by us- 
ing a stochastic absorption model. Recent detailed nu- 
merical studies have been performed by transfer matrix 
methods [20,21] to study our approach, and the numer- 
ical results support our model. Here, we first study 
the coherently absorbing disordered media modeled by 
a adding complex potential. Then, we present a detailed 
derivation and analysis of the Langevin equation in case 
of stochastic absorption for the modeling of absorption 
by "fake" , or "side" , channels obviating the need for a 
complex potential. We show that stochastic absorption 
model gives potentially better physical results than the 
model by a complex potential for transport in ID absorb- 
ing disordered media in the regime of strong absorption. 
However, these two models give the same result for weak 
disorder and the weakly absorbing regime. 



II. COHERENT ABSORPTION MODELED BY 
COMPLEX POTENTIL 

There has been recent interest in the role of absorption 
in localization, and on the different length scales in the 
problem in the presence of absorption. For the bosonic 
case, like light etc., a coherent state (e.g. a laser beam) 
is an eigenstate of the annihilation operator. Removal of 
a photon (absorption) does not destroy the phase coher- 
ence. Fermions cannot be annihilated in the context we 
are considering here. For the Fermionic case, the phys- 
ical picturization of the absorption will be some kind of 
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inelastic process (like scattering by phonons) , where elec- 
trons lose their partial temporal phase memory. This 
type of absorption is called stochastic absorption. 

Recent theoretical studies [22-24] have shown that the 
absorption in the case of light waves does not give any 
cut-off length scale for the localization problem. If a sys- 
tem is in the localized state, absorption will not kill the 
localization to make the system again diffusive. A sharp 
mobility edge exists even in the presence of significant 
absorption in 3D. 

Adding a constant imaginary part with the proper 
sign to the real potential of the Maxwell/Schrodinger 
wave equation can model the linear coherent ab- 
sorption/amplification. Both the Maxwell and the 
Schrodinger equations can be transformed to the form of 
the Hclmholtz equation. A general form of a Helmholtz 
equation can be written as: 



^ + k''[l + {r]{x)+ir]a)]u 



= 



(1) 



For the electronic case, consider V{x) + iVa as the po- 
tential, k as the wave vector, and Ji/2m = 1. For the 
optical case, consider eo -I- zEq as the constant dielectric 
background and e{x) as the randomly spatially fluctuat- 
ing part of the dielectric constant, k = 2n/X, where A = 
wavelength in the average medium (eg). Then, we can 
define rj{x)——V{x)/k'^, r]a=—Va/k^ for the Schrodinger 

equation, and r}{x)= "'"f^' ^"^^ = for 

the Maxwell equation. 

The model we are considering here is for a one-channel 
scalar wave where the polarization aspects have been ig- 
nored. This holds well for a single-mode polarization 
maintaining optical fiber. For electronic case, we con- 
sider only a single channel. 

As pointed out in Ref [19], modeling absorption by a 
complex potential will always give a concomitant refiec- 
tion. Absorption is not possible without reflection and 
there is a competition between the absorption and the 
reflection. For very high absorption, the system may try 
to act as a perfect reflector. We extend and explore this 
idea for coherently absorbing 1-D disordered media. The 
Langevin equation for the complex amplitude reflection 
coefficient R{L) can be derived by the invariant imbed- 
ding method [25] from the Helmholtz equation Eq.(l), 

= 2ikR{L) + z J(l + 77(L) + ir^a)il + R{L)f , 

(2) 

with the initial condition R(L)=0 for L=0. 

Now, taking R{L) ~ -y/re*^, the Langevin equation 
Eq.(2) reduces to two coupled differential equations in 
r and 6. For a Gaussian white noise potential, using a 
stochastic Liouville equation for the evolution of proba- 
bility density and then integrating out the stochastic part 



by Novikov's theorem, we get the Fokker-Planck equa- 
tion: 



dP{r, e) 
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(3) 



where we have introduced the dimensionless length I = 
and D = ^ = ^/^a, and ^ = {\qk'^)~^ is the localization 
length and = ^^"^ absorbing length scale in the 

problem, and q is tlie strength of the delta correlation of 
the random potential. 

The Fokker-Planck equation Eq.(3) can be solved an- 
alytically for the asymptotic limit of large lengths in the 
random phase approximation (i.e., for the weak disorder 
case). Then, the statistics P^{r,l) have a steady state 
distribution for Z ^ oo. It can be obtained by setting 
dP{r,l)/dl = and solving the resulting equation ana- 
lytically. We get [4] the following steady state distribu- 
tion: 



lim P^(r,Z) = 



\D\exp{\D\)exp{-j^) 





for r < 1 
for r > 1, 



(4) 



For D = and I ^ oo, equation (3) reduces to a 
half delta function peaked around r = 1. The Fokker- 
Planck equation, Eq.(3), is difficult to solve analytically 
beyond the random phase approximation. Wc solve the 
full Fokker-Planck equation (3) rnimcrically to calculate 
the statistics of the reflection coefficient and its phase in 
a regime of strong disorder with strong absorption. 

To sec the actual evaluation of the probability as 
a function of reflection coefficent for a fixed disorder 
strength and length, we evolved the full Fokker-Planck 
equation without approximation. The Fokker-Planck 
equation is solved as an initial value problem, i.e. P(r, I = 
0) =S{r = 0). 

In Fig.l we plot the results of our numerical simula- 
tions. The P(r) distributions show (top plot) that for 
the strongly absorbing disordered media, the probability 
initially peaks near r = and the peak slowly moves to- 
wards r = 1 for larger sample lengths, indicating stronger 
reflection for stronger absorption. In Fig.l, the bottom 
plot shows that the steady state phase distribution is a 
delta function peaking symmetrically about 6 = tt. 

These plots, show the system behaving as a perfect re- 
flector for higher lengths, with strong disorder and strong 
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FIG. 1. Probability distribution P{r) and P{0) separately 
against the sample length I. Parameters are fixed strong dis- 
order strength 2kF^ = -001 and fixed strong absorbing par- 
rameter \D\ = 10. 

absorbing parameter Therefore, this coherent ab- 

sorbing model is restricted for the weakly absorbing pa- 
rameter regime. 



III. STOCHASTIC ABSORPTION 

In the previous section, we have shown and discussed 
that the absorption is not possible without reflection in 
a calculation where absorption is modeled by adding a 
constant imaginary potential to the real potential. In the 
Langevin equation for R{L), derived from these types of 
models, the wave always gets a reflected part along with 
the absorption. There is, however, another way of deriv- 
ing a Langevin equation for the reflection amplitude R{L) 
such that the absorption does not have a concomitant re- 
flected part. This approach is motivated by the work of 
references [26,27] where some purely absorptive "fake" 
(side) channels arc added to the purely clastic scattering 
channels of interest. A particle that once enters to the 
absorbing channel, never comes back and it is physically 
lost. We derive the Langevin equation for R{L) follow- 
ing the approach in Ref. [26]. This Langevin equation 
has some formal differences from Eq.(2). However, 




Reservoir 

FIG. 2. Modeling "absorption" by fake channels; Channels 
1 and 2 are coupled through the current leads to two "fake" 
channels 3 and 4 which connect to a thermal reservoir. 

it does not differ for the weak disorder case. 

Now we proceed to derive Langevin equation for R{L) 
using the "fake" channel approach. The main idea is to 
simulate absorption by enlarging the S-matrix to include 
some "fake" , or "side" channels that remove some prob- 
ability flux out of reckoning. Let us assume that the 
scatterer has a general scattering channel and an absorb- 
ing channel ending with a phase randomizer which acts 
as a blackbody. The scattering matrix can be taken as 
discussed in reference [26]. Consider scattering channels 
1 and 2 connected through current leads to two quan- 
tum "fake" channels 3 and 4 that carry electrons to the 
reservoir, ( with chemical potential /x) as shown in Fig. 2 
(for only one scatterer). This is a phenomenological way 
of modeling absorption [26,27]. Electrons entering into 
the channels 3 and 4 arc absorbed regardless of their 
phase and energy. The absorption is proportional to the 
strength of the "coupling parameter" e. A symmetric 
scattering matrix S for such a system can be written as : 



( aR qT 

aT aR 







f3 \ 








aR* 


-aT* 


aT* 


-aR* J 



(5) 



where R and T are the reflection and the transmission 
amplitudes for the single scatterer , a = \/\ — t is the 
absorption coefScient, and /3 = \fk. 

Now, the full scattering matrix is unitary for all pos- 
itive real e < 1. But the sub-matrix directly connect- 
ing only the channels 1 and 2 is not unitary. We will 
now explore this fact in our derivation. We observe from 
the above sub-matrix that for every scattering involving 
channels 1 and 2 only the reflection and the transmission 
amplitudes get multiplied by a (the absorption parame- 
ter). Now keeping this fact in mind, we will derive the 
Langevin equation for the reflection amplitude R{V) for 
N scatterers (of length L) each characterized by a ran- 
dom S-matrix of this type, with the same a value. 
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FIG. 3. A schematic picture of a scattcrcr of length L hav- 
ing reflection amplitude R[L) while an added length dL makes 
the length of the scatterer L + dL. The effective potential of 
the length dL is shown by an effective summary deltar function 
potential. 

Consider N random scatterers for a ID sample of 

length L with (= TVa, a = unit-spacing) the reflection 
amplitude R{L), and + 1 scatterers of the sample 
length L + AL, with the reflection amplitude R{L + AL). 
That is, let us start with the N scatterers, and add one 
more scatterer to the right to make up the N+l scat- 
terers. Now, we want to see the relation between R{L) 
and R{L + AL). That is, let there be a delta-function 
potential scatterer between L and L + AL, positioned at 
the point L + AL/2, that can be considered as the eff'ec- 
tive scatterer due to the extra added scatterers for length 
AL. For kAL <C 1, we can treat the extra added scat- 
terer as an effective delta potential vq {L)6{x — L — AL /2) 
with Vo{L) = V(L)AL.{Wc consider the continuum limit, 
a ^ 0, ^ oo, and Na = L fixed ). 

Now, for a plane-wave scattering problem for a delta 
hmction potential of strength vq, which is at .t = and 
has complex reflection and transmission amplitudes ro 
and to respectively, we have from the continuity condi- 
tion for the wave function and discontinuity condition 
for the derivative of the wave function (which one gets 
by integrating the Schrodinger equation across the delta 
function): 



ro 



i2k — vo 



and to = ro + 1- 



Considering Wq = v{L)AL the smallness parameter, one 
gets expressions up to first order in AL for ro and to- 

ro = ^ and ^0 = 1 + ^, 

where we have taken jlm = 1. 

Now to introduce absorption we will write: 

r— > ra = ro(l — e)^^^, 
t^ ta = to(l - e)^/^- 



This means that for every scattering the refiection and 
transmission amplitudes get modulated by a factor of a. 

Now consider a plane wave incident on the right side of 
the sample of length L + AL- Summing all the processes 
of direct and multiple refiections and transmissions, on 
the right side of the sample of length L, with the effective 
delta potential at L + AL/2, we get, 

R{L + AL) = re^^^^ 

_|_ gifeAL/2^gjfcAL/2^^2;^'jgi/cAL/2^gjfcAL/2 

+ (6) 

Summing the above geometric scries, substituting the 
values of r and t, and taking the continuum limit for 
L , one gets from Eq.(6), the Langevin equation: 



dR{L) 
dL 



-aR{L) + 2ikR{L) + i-(ry(L))(l + R{L)f , 



(7) 



with the initial condition R{L) = for L = 0, and a is 
the absorbing parameter. This is not quite the same as 
the Eq.(2) obtained by introducing an imaginary poten- 
tial. It turns out, however, that this Langevin equation 
gives the same results in the regime of weak disorder but 
differs qualitatively in the regime of strongly absorbing 
disordered media. 

Similarly, following the same steps as we discussed pre- 
viously, from the Langevin Eq.(6) we derive the Fokker- 
Planck equation: 



dP{r, 9) 

di 



sm6'-^(r^/^(l - r)) ^ 



dr 
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+ \{r"'+r"')§-^{cosm 



P{r, 6) 



- 2k^ 



dP{r, ( 



+ D 



d{Pr) 



(8) 



86 dr 

where parameters I and 2fc^ are same as defined in Eq.(3), 
and £> = 

In the weak disorder limit, when the phase part is in- 
tegrated out by the RPA, Eq. (8) and Eq.(3) give the 
same Fokker-Planck equation in r. Hence, the steady 
state P{r) solution for the both are same, i.e. Eq.(4). 
We will consider here only the strong disorder case for 
the numerical calculations of the Fokker-Planck equation 
(8), which goes beyond the RPA to compare it with the 
solution of Eq.(3). 

Fig. 4(a) is the plot of P{r) for the strong disorder 
{2kF^ = .001) and strongly absorbing {D = 10) limit, 
with different lengths of the sample. P{r) distributions 
are peaked around r — 0. The wave is absorbed in the 
medium before it gets reflected as we have modeled ab- 
sorption without reflections and we arc considering the 
case of strong absorption, and hence the probability of 
absorption is more than that of reflection. 



4 



<D 
T3 



CO 
O 



(a) 



2k5=.001 
D=10 



1=1,5 



as a perfect reflector. Therefore, tlic Langcnvin equation 
derived from stochastic absorption method gives a poten- 
tially better physical model for coherent absorption. Our 
model proposed here has already been verified and sup- 
ported by extensive numerical simulation using a transfer 
matrix method [20,21]. 

I gratefully acknowledge N. Kumar for many stimulat- 
ing discussions. 



.01 .5 1 

Reflection coefficient r 




Phase 



FIG. 4. Probability distributions (a) P{r) and (b) P{e) 
against the sample length I with fixed strong disorder strength 
2kF^ = .001 and fixed strong stochastic absorption parameter 
D = 10. 

Fig. 4(b) shows the phase distribution P{9), which is 
a double peaked symmetric distribution and differs from 
the model of absorption by complex potential, which is 
a delta function distribution at ^ = tt as shown in Fig.l 
(bottom). 

IV. DISCUSSION AND CONCLUSION 

In conclusion, we have shown a comparison between 
the modeling of absorption in disordered medium (i) by 
a complex potential and (ii) by fake channel S matrix 
approach. We derived a new Langcvin equation using 
the latter method. We calculated the statistics of the re- 
flection coefficient (r) and its associated phase (6) for a 
wave reflected from absorbing disordered media, for dif- 
ferent disorder strengths and lengths of the sample, for 
the both models. For the case of weak disorder with 
weak absorption, i.e. within the RPA, the FP equation 
for the case of stochastic absorption is the same as that 
of the absorption modeled by a complex potential. For 
the regime of strong disorder with strong absorption, the 
stochastic absorption model behaves as a perfect absorber 
opposite to the coherent absorbtion case, which behaves 
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